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^ ■ Abstract 

V ' Anomalous production of low-energy photons from the galactic center have fueled 

speculations on the nature and properties of dark matter particles. In particular, it 
r~| I has been proposed that light scalars may be responsible for the bulk of the matter 

density of the universe, and that they couple to ordinary matter through a light 
^Sj . spin-1 boson. If this is the case, then such particles may be produced in the quasi- 

\ elastic low-energy scattering of electrons off protons. We present a proposal for an 

' experiment to search for this process and assess its viability. 

o " 
in 

o ; 1 Introduction 

!>■ 
O 



X 



The INTEGRAL data provide evidence for an anomalous production of low-energy gammas 
from the galactic center [1]. The excess is largest in the region at and below = 511 keV 
leading to the conclusion that positron-electron annihilation is responsible. The question 
becomes: where are the extra positrons coming from? One answer is that they are produced 
in the annihilation processes of dark- matter particles, which we will write as XX* ~^ e^e~ , 
where x denotes the dark-matter (DM) particle. 

A phenomenological model has been developed which posits a light dark-matter particle x 
and a new spin-1 gauge boson U which mediates interactions between x and ordinary matter 
particles, such as electrons [2]. The couplings of the ?7-boson to standard model fermions 
is not specified by any theoretical first principles; rather, they are left as parameters to be 
constrained by the known dark-matter abundance and processes studied in the laboratory. 
At a minimum, there must be a non-zero coupling to the dark matter particles as well as to 
electrons, a fact which we exploit in our proposal to observe U bosons in the laboratory. 

The experiment we propose makes use of a low-energy electron beam (-Ebeam ~ 40 MeV) 
and a fixed hydrogen gas-jet target. We describe two versions: the first is rather simple and 
requires a minimum of resources, and the second is more elaborate allowing for a greater 
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sensitivity. We attempt to compare the sensitivity and potential information gained from 
our proposed experiment to other suggested avenues of research, such as the study of e~^e~ 
colhsion data [3]. 

This paper is organized as follows. We begin with a brief review of the evidence for the 
low-energy gamma-ray excess from the galactic center, followed by a resume of the light dark 
matter models recently discussed in the literature. Next we explain the process of interest 
in general terms, followed by a detailed description of the two experimental set-ups to study 
this process. We make a comparison with other proposals followed by a summary of our 
work. 

2 Gamma- Ray Excess 

In 1972, Johnson et al. [4] used a balloon-hoisted Nal scintillation telescope to detect a 
511 keV line emanating from the center of the galaxy. After accounting for 511 keV radiation 
resulting from cosmic rays and positron annihilation in the upper atmosphere, they concluded 
that this line was significantly stronger than the cosmic background radiation from other 
directions at similar energies. Leventhal et al. [5] revisited this phenomenon in 1978 and 
identified the source of the radiation as positron annihilation; specifically, they concluded 
that positrons in the galactic center form both para-positronium, which annihilates to give 
the 511 keV line, and ortho-positronium which contributes a spectrum of excess low-energy 
radiation at and below 511 keV. Further experiments over the past 30 years, and most 
recently the INTEGRAL satellite, have refined the value of the observed photon flux and 
identified the majority of the radiation as coming from the galactic bulge [1,6-8]. Currently, 
the explanation of annihilating positronium is well-accepted as the source of the radiation, 
but the source of the positrons themselves remains unclear. Several possibilities involving 
relatively well-known astrophysical phenomena have been put forth, including radioactive 
nuclei from supernovae [9], gamma-ray bursts [10], pulsars [11], black holes [12], and cosmic 
rays [13], but these models have difficulty accounting for the morphology and high intensity 
of the photon flux except under rather restrictive assumptions. 

To resolve these problems, several more exotic explanations from particle physics have 
been proposed. A partial list compiled by Sizun et al. [14] includes Q-balls, relic particles, 
decaying axinos, primordial black holes, color superconducting dark matter, superconduct- 
ing cosmic strings, dark energy stars, moduli decays, and annihilating light dark matter 
particles. The latter explanation, in which dark matter particles annihilate in the galactic 
bulge to form e"'"e~ pairs, has been able to account for both the morphology and intensity of 
the line as long as the dark matter particles are light (< 100 MeV) [15-18]. Radiative pro- 
cesses play an important role in comparisons with gamma-ray data, and provide additional 
constrants on the model [14,19]. In the past four years, the light dark matter model has 
been refined to account for the increasingly accurate data from INTEGRAL; the upper limit 
of the particles' mass has been placed anywhere from 3 to 20 MeV [14,20,21,23] and the 
theoretical aspects of this model have been studied considerably [17,24,25]. In particular, the 
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flux intensity measurements have allowed a fairly accurate determination of the cross-section 
of the annihilation reaction [25,27]. Several authors have noted that increasingly accurate 
measurements of the morphology of the 511 keV line would allow the mass of the proposed 
particle to be narrowed down even further [20-22]. 

The annihilating dark matter scenario has attracted much attention because it is the eas- 
iest "exotic" explanation to test experimentally. For instance, Ref. [28] focused on detecting 
a line from the galactic center that would result from the dark matter particles annihilating 
directly to photons, and Ref. [3] proposed searches for the gauge boson involved in annihi- 
lation processes at high-energy e^e~ colliders. We show that the light dark matter model 
can be tested cleanly and simply in the laboratory, using techniques and methods that are 
readily available today. Even taking into account the range of possible dark matter particle 
mass, the kinematical distinctions between a dark matter production event and background 
events are clear enough that a signal could be identified with a high degree of confidence. 
Moreover, if dark matter were to be detected in such an experiment, the measurement of the 
cross-section would allow an independent check of the parameters in the light dark matter 
model, and by extension, provide significant information on the annihilations that take place 
at the center of the galaxy. 

3 Resume of the Light Dark Matter Model 

If the positron excess is explained by XX* ~^ e~^e~, then clearly there is an effective e'^e~xx* 
vertex. Calculations by Boehm, Fayet and others have linked the strength of this effective 
vertex to the observed dark matter abundance, and the dark matter particle mass, M^. The 
e'^e~xx* interaction would be modeled by one or more intermediate particles. According to 
a recent paper [25] , one needs both a heavy fermion called the and a light neutral vector 
boson called the U to explain both the primordial abundance of dark matter and the current 
rate of positron annihilation in the galactic center. In this model, the dark matter particles 
are neutral scalars with masses on the order of 1-10 MeV. There are other models in which 
X is a fermion [26], but for concreteness, we will take x to be a scalar. The f/-boson is also 
light, with M[/ > 2M^ and M[/ < 100 MeV, and it decays mainly to xX*- The mass of the 
fermion is at least several hundred GeV and plays no role in our process, so we do not 
consider it further. Feynman diagrams for the dark matter annihilation are given in Fig. [TJ 

The experiment that we propose does not depend crucially on these assumptions. Rather, 
this phenomenological model serves as a guide for gauging the sensitivity of our experiment. 
In order to estimate rates for our signal process, we must specify the coupling constants. 
We follow a recent paper by Fayet [29] , which links this model to several experimental data, 
as well as a paper by Ascasibar, et al. [25], which links the model to various astrophysical 
data. According to Ascasibar et al., the U boson plays the dominant role in fixing the 
primordial abundance of dark matter, while the fermion controls the present day rate 
of annihilations, and hence the rate of positron production. Fayet uses the dark matter 
abundance to constrain the ?7-boson couplings and does not hypothesize a heavy fermion. 
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Figure 1: dark matter annihilation into an e'^e pair, through U -boson exchange in the 
s-channel, and F^-fermion exchange in the t-channel 



He shows how particle physics data constrain the U -boson couphngs as a function of its mass. 
Despite these differences, both papers report similar constraints on the coupling constants 
of the U boson to the dark matter particles, C^, and to electrons, fe- Other publications 
placing similar constraints include [2,24,30,31]. 

If the dark matter abundance is used to constrain the f/-boson couplings, then 

'^^^^'"^^^ MeV) V^^"-^ 

where -Bann is the fraction of all XX* annihilations which result in an e~^e~ final state (see 
Eq. (57) in Ref. [29]). We will assume that -Bann — although a lower value is possible if 
there is a significant coupling of the ?7-boson to neutrinos. If x were a spin-1/2 fermion, then 
this formula would be modified by a factor 0{1). Since our aim is to assess the viability 
of an experiment, however, we will neglect any such factors in this paper. Fig. [2] (TOP) 
depicts this constraint as a function of for three choices of Mf/. The single dot shows our 
default choice of model on which our rate estimates in Sec. Hlare based. Fig. [2] (BOTTOM) 
shows the constraint from Eq. ([1]) as a function of Mu, for three choices of M^. It also 
shows the upper bounds on |C^/e| derived by Fayet from three lepton-based experimental 
measurements, namely, the measurement of {g — 2)e [32], the measurement of {g — 2)^ [33], 
and the u — e cross-section measurement [34]. These constraints apply to the vector coupling 
of the f/-boson to the electron; more stringent constrains apply to the axial coupling. We 
will set the axial coupling to zero in our calculations. The constraint from u — e scattering 
assumes that f,^ = fe and will be weaker if <C fe- The constraint from {g — 2)^ applies to 
our process only for = /g. Fayet also derives constrains from ip and T decays, as well as 
atomic parity violation experiments, but we will assume that the t/-boson couples to leptons 
only. 
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Figure 2: TOP: constraints on C^fe as a function of M^, for three values of Mu. BOT- 
TOM: heavy red lines show constraints on C^fe as a function of Mu, for three values of M^. 
The thinner dashed blue lines show constraints coming from the measurements of {g — 2)^, 
{g — 2)e and a{ve) (see text), from Ref. [29], and we have set = 1. The hook-shaped green 
curve indicates the sensitivity of the first experimental design, discussed in Sec. \5.1[ The 
dots indicate our default choice for masses and couplings. 



6 



Low-Energy ep Scattering 




Figure 3: example Feynman diagrams for the signal process. We consider only the U-hoson 
diagrams in our estimates; the F^-fermion diagrams can he neglected. 



4 Low-Energy Quasi-Elastic e;?- Scattering 

We wish to exploit the different kinematic characteristics of the signal process 

e~p e~pU^*^ ^~PXX* 

and of simple elastic scattering 

e~p e~p 

at low energy. The kinematics of the scattered electron and proton are highly constrained 
for elastic scattering, and much less so for the signal process. Several kinematic distinctions 
can be made, which allows an efficient and effective discrimination of the two processes on 
the basis of kinematics alone. 

Feynman diagrams for the signal process are given in Fig. [31 It is important to note that 
the same set of vertices appear in both our signal process and the dark matter annihilation 
process (Fig. [T]). Consequently, our signal process must occur if this model is correct, and 
the signal cross-section can be related directly to the dark-matter annihilation rates in the 
early universe (which determines the dark matter relic density) and today (which determines 
the strength of the 511 keV gamma-ray line from the galactic center). 

Since the dark matter particles are light, by hypothesis, we choose to employ a low 
electron beam energy which essentially closes off all inelastic processes, leaving only elastic 
scattering as a significant background process. For example, we may take an electron beam of 
energy -Ebeam ~ 40 MeV and a fixed hydrogen target, for which y/s ~ Mp+Ebeum < Mp+Mt^, 
since -Ebeam ^ Mp. There can be no inelastic background from ordinary strong- interaction 
processes, and if a pure hydrogen target is used, there would be no nuclear excitations, either. 

A very small background will come from the process e~p e~ph'u mediated by an off- 
shell Z boson. The Feynman diagram is the same as the ?7-boson exchange diagram (Fig. [3]), 
with the t/-boson replaced by a virtual Z-boson and the scalar particles XX* replaced by uu. 
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This background process is completely negligible - about nine orders of magnitude smaller 
than the hypothetical signal - since Mu -C Mz- 

Higher-order QED radiative processes e~p ^~Pl{l) also pose a background. They 
occur at a smaller rate than elastic scattering. The final-state photons tend to emerge 
along the directions of the incoming and outgoing electrons, and hence do not change their 
directions [35], although a small component of "wide-angle" Bremsstrahlung can occur and 
has been observed in muon-scattering experiments [36]. Calculations for experiments at 
HERA and JLAB show extremely small rates outside a cone of 0.1 rad [37]. In the case that 
the radiated photons emerge along the scattered electron direction, they can be detected 
together with the scattered electron. If they are emitted along the beam direction, thereby 
reducing the effective ^/s for the ep interaction, the correlation of the electron and proton 
directions changes very little. This is a special feature of our kinematical situation, in which 
Mp ^ -Ebeam ^ ^e- Additionally, QED radiation peaks at low photon energies, which means 
that the deviation of the scattered electron energy from the lowest-order results tends to be 
small, with no sharp kinematic features. The signal process peaks at large deviations, and 
there are distinctive kinematic features due to the phase space needed to produce two dark- 
matter particles, and/or an on-shell t/-boson. Finally, the QED radiation pattern varies very 
little with the beam energy, while the signal process will have a strong dependence. Hence, 
discrimination between the signal process and QED radiative processes will be possible if 
the signal rate is not too small. 

The concept for the experiment is the following: collide a well-defined electron beam onto 
a hydrogen target, and observe the scattering angle and energy of the outgoing electrons 
and protons. For ordinary elastic scattering, the energy of the scattered electron emerging 
at a given angle is constrained to a unique value. So one would look for events with a 
significantly lower energy as evidence for the production of a final state that is lighter than a 
single pion. Confirmation for an anomalous final state would come from the measurement of 
the scattering angle and energy of the outgoing proton. The kinematic distributions for the 
final-state electron and proton would allow, in principle, a confirmation of the production 
of light invisible scalars of a particular mass, as opposed to the production of a pair of 
neutrinos. In the absence of a signal, stringent limits could be placed on the masses and 
effective couplings of dark matter to electrons. 

The keys to identifying signal events are: 

1. For a given scattering angle, the scattered electron will have a much lower energy than 
in elastic scattering. 

2. The outgoing proton will be relatively slow; for elastic scattering the proton is energetic. 

3. For a given electron scattering angle, the scattering angle of the proton will vary over 
a wide range. In elastic scattering, the proton emerges at a unique angle. 

4. The electron and proton can be acoplanar for the signal event, due to the momentum 



8 



Low-Energy ep Scattering 



taken by the XX* pair; for elastic scattering the electron and proton are strictly back- 
to-back in the transverse plane, even if there is final-state radiation. 

5. The signal will increase rapidly as -Ebeam is increased from threshold, while elastic 
scattering decreases as There will be a minimum beam energy corresponding 
to the threshold for producing two dark-matter particles, below which there is no signal. 

6. The signal cross-section is less sharply peaked toward small electron scattering angles. 

7. In principle, the U boson may decay to an electron-positron pair rather than the 
invisible xX* state. The final topology would contain four charged tracks. If the U- 
boson is on mass shell, then the mass distribution of the 'extra' e~^e~ pair would give 
a peak at the f/-boson mass. If it is off mass shell, then the mass distribution will be 
less sharply peaked toward 2Me than that due to photon conversions. 

These very distinctive features allow an efficient selection of events with very little back- 
ground. As described later in this paper, a simple apparatus should allow signal events to be 
identified at the rate of one event in a ten billion elastic scatters, and a more sophisticated 
apparatus should achieve a sensitivity better than 10~^^ of the elastic scattering rate. Such 
experiments should easily cover the range of possible models of this type, leading either to 
the discovery of the [/-boson or the exclusion of this and similar models. 

While we have based our calculations on a fairly specific model and final state, our 
argument is not dependent on this model in all of its details. As already noted, the x particle 
might be a fermion instead of a scalar. Furthermore, the experiment we propose could be 
viewed as the direct production of ?7-bosons, regardless of how they decay or whether they 
play any role in dark matter phenomena. In this case, the rates would not depend at all 
on the coupling C^, and the final states might be dominated by other light particles. If a 
signal for an anomlous invisible final state were observed, then follow-on experiments would 
be needed to confirm the connection with dark matter (see Ref. [38]). In order to provide a 
clear and consistent framework for our discussion, however, we will follow fairly closely the 
light dark matter model described above. 

We proceed now to a discussion of the rates. 

4.1 Kinematics 

We remind the reader of the basic kinematics for elastic scattering in order to frame our 
discussion and define our notation. The target proton is effectively at rest in the lab- 
oratory frame with four-momentum = (M, 0,0,0), and the incoming beam electron 
is highly relativistic with four-momentum = {E,0,0,p). The outgoing electron has 
p'^ = [E' ,p' sm6,0,p' cos6) which defines the electron scattering angle 6. The outgoing 
proton has P'^ and the four- vector of both dark-matter particles we will write as W^. Hence, 
Wfj, = p^ + — p'^ — P'^, and = for elastic scattering. For the signal process, 
mm{Wf,Wf') =4:M^. 
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Following decades-old convention, we define the four- momentum transferred = p'^—Pn 
and Q2 = _g2 ^ very good approximation, = AEE' sin^ ^/2. The energy transfered in 
the lab frame is u = E — E' . The kinematic condition for elastic scattering, P'^ = = M^, 
implies ~ 2Mu, after neglecting a small term proportional to m^. 



4.2 Elastic Cross-Section 

The elastic electron-proton scattering cross-section, for relativistic electrons is 
da a'^n^c^ ( E'' 



dVL 4E2sin^(^/2) 
where r = Q'^/AM'^ and 



i + r 



are functions of the proton electric form factor Ge and magnetic form factor Gm. At our 
low beam energies, r ^ 1, so we can neglect the magnetic form factor. G^; is approximated 
by the standard dipole fit 

where /? = 710 MeV. Under these conditions, Eq. ([2]) simplifies to 

da a^h^c^ (E'\ ( 



dn AE^sin\e/2) \ E ) + Q\ 
E' and 9 are related by 



cos\e/2). 



1 + {E/M){1 - cost 



The contribution of the form factor to the cross-section is on the order of a few percent at 
low beam energies; this correction is necessary to get an accurate estimate of the elastic 
background, but is not needed for an order-of-magnitude estimate of the signal process. 



4.3 Signal Cross-Section 

We investigated the signal process in two ways. First, we performed a semi- analytic calcula- 
tion in which the matrix element was assumed to be independent of all momenta. Then we 
employed CompHEP [39] and implemented the ^/-exchange process, an example of which is 
shown in the left-hand Feynman diagram of Fig. [31 
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4.3.1 Semi-Analytical Cross-section 

To get a rough idea of the signal kinematics, we calculated cPa / dVLdE' assuming that the 
matrix element of the signal process is independent of all momenta. The phase-space calcu- 
lation for a 2 — > 4 process can be rather involved, especially when several of the final-state 
particles are massive, and finding an analytic expression for d'^a/dQ dE' in terms of the 
four- vectors of all outgoing particles is not possible in general. We simplified the calculation 
by writing VK^ for the sum of the two four-vectors of the scalar particles and expressing the 
cross-section in terms of the combined invariant mass and the combined spatial velocity 
W, treating these two quantities as separate variables. In the calculations that follow, all 
energies and momenta pertain to the center-of-mass frame. 
We begin with the Golden Rule for 2-^3 scattering: 



da = \M\' 



d^P' \ / d^W 



d^p' \ 

{2Txf2E' ) \{2Tif2E' ] \{2Txf2E 



x(27r)V(j9 +P -p' -P' -W) 

where the relevant four- vectors are defined in Section 14. 1[ Since we are only concerned with 
the shape of the phase space distribution, we will drop all numerical constants for the rest 
of this calculation. Making the standard approximation of a massless electron, we integrate 
out the proton momentum and the angular part of W, using the delta function to obtain 

— * 

limits on \W\: 



_ EE' ± (Etot - E')^F^ - 4iy2^tot(^tot - 2E') 
2£'tot(-E'tot — 2E') 

where Etot is the total energy, E' is the outgoing electron energy, and F = E^^^ + — 
2EtotE' - M^. At this point 

rfV r\^\+ \W\ 



dQ.dE' J\w\. Iy^2 _^ 1^12 



oc ' d\W\ = ^W'^^\W\\-^W^^\W\^_. (4) 



Requiring that |Vr|± be real, i.e., that the term under the radical in ([3]) be non-negative, 
gives the upper bound 



{W)^ = M^ + El, - 2E,,,E' - 2M^El, - 2E,,,E' 

with the lower bound {W'^)^ = 4:M^. We use these bounds to integrate the right-hand side 
of Eq. (j4]) over W^. It turns out that this integrand can be very well approximated by a 
square root function: 



f{W^) = + \W\l - \JW^ + \W\l ^ KJ{W^)+ - W^. 
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Figure 4: Feynman rules for the U and U -e vertices 



The normalization constant K is chosen so that the square root function matches f{W^) at 
the endpoint = (W^)-; that is, 



(The other endpoint Pt^^ = is aheady taken care of since both the square root and 

fiyV^) vanish there.) Performing the integration over VT^, we now have 

2^ AW^-\. flf^M2\ \ 



d^a f(w')+ f{{W 



dn dE' 



oc / , ■^^^ J{W^)+-W^d{W 



oc f{{W^).){{W^)^-AM. 



X' 



AMI + 1^1+ " V^^^ + ((^ )+ - 4M^) 

where \W\± are evaluated at = 4M^. 

A plot of this function, boosted from the CM frame to the lab frame, gives the broad 
phase-space curves in Fig. [71 We note that this calculation is model-independent in the sense 
that it makes no reference to the nature of the interaction. 



4.3.2 CompHEP Calculations 

We implemented [/-boson exchange between electrons and dark matter particles. An example 
Feynman diagram is shown on the left side of Fig. [31 The proton was represented as a massive 
fermion with no internal structure, which, for the very low beam energies we propose, is a 
very good approximation. The mass of the dark matter particle was set to M^ = 2 MeV, 
and the mass of the ?7-boson was set to Mu = 10 MeV. Hence, the ?7-boson was produced 
on mass shell. The coupling of the [/-boson is assumed to be purely vectorial - see Fig. [Hfor 
the Feynman rules for the U-x and U-e vertices. We imposed the relic abundance constraint, 
Eq. dH). Taking C^^ = I and B^^^ = 1, this gives fe ^ 2.3 x 10"^ 

For the purposes of this paper, we set the beam energy to -Ebcam = 40 MeV, and we 
considered a narrow range of electron scattering angle 89.4° < 9^ < 90.6°, which follows 
from the conceptual design presented in Section 15.11 
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The mass of the ?7-boson is not well constrained. If the product of coupling constants 
C)(^fe is held fixed, then the signal cross-section drops rapidly as a function of Mu, as depicted 
by the solid curve in Fig. [5] (TOP). The presence of the threshold at Mu ^ 40 MeV is evident. 
If the Mu boson were heavier, then one would run this experiment at a higher beam energy. 
Backgrounds will not increase so long as -Eboam < 135 MeV, and in fact the signal cross- 
section will increase, up to a certain point. According to Eq. ([1]), the constraints C^fe should 
be adjusted as a function of Mu- The result for the cross-section is shown as the dashed 
curve, which is much flatter than the solid curve, showing that the rate of this process is 
indeed tied directly to the rate of dark matter annihilation. Fig. \5\ (BOTTOM) shows the 
variation of the signal cross-section as a function of beam energy, for three choices of Mu- 
A typical threshold behavior is evident. 

Since the signal process is inelastic, the electron is scattered at larger angles than for 
the background process - see Fig. [6] (TOP). Clearly the ratio of signal to background is a 
strong function of the scattering angle, and so we have chosen 6 = 90°, as an example. An 
optimized choice of 6 would depend on Mu as well as the details of the apparatus, which is 
beyond the scope of this paper. 

The wide range of proton scattering angles, 6p, is also a distinctive feature of the signal 
process. Fig. E] (BOTTOM) shows da/d6p, with the electron scattering angle 6e constrained 
as above. The dashed lines indicate the narrow range of 6p expected for elastic scattering. 
This range does not change even when the incoming electron emits an energetic photon 
during the scattering process. 

An important feature of the signal process is the reduced energy of the outgoing electron. 
For -Ebeam = 40 McV and 6 = 90°, perfect elastic scattering gives E' = 38.4 MeV. For the 
signal process, E' follows a much broader distribution with a peak at rather low values, as 
shown in Fig. [71 It is an interesting feature of quantum field theory that a massless spin- 
1 boson, such as the photon, is emitted with an energy that peaks toward the minimum 
possible value, while a massive spin-1 boson, such as the [/-boson considered here, is emitted 
with a momentum that peaks toward the maximum possible value. Thus the CompHEP 
calculation, based on the matrix element for [/-boson exchange, shows a scattered electron 
energy distribution which peaks toward small values. The phase-space calculation, described 
in Section 14.3.11 sets the matrix element to a constant, resulting in a broader distribution 
for the scattered electron energy as shown in Fig. [3 

On the basis of these CompHEP calculations, and for 89.4° < 6^ < 90.6°, the accepted 
signal cross-section is 2.13 x 10~^ pb, and the cross-section for elastic scattering is 8.0 x 
10^ pb, which gives a ratio of cross-sections of 2.7 x 10^^°. We checked these results using 
CalcHEP [40]. In the next section we describe experiments which should be able to probe 
this range successfully, and either establish the existence of this signal process or rule it out 
definitively. 



An Experiment to Search for Light Dark Matter 



13 



1 



o 

• I— ( 
^— > 

o 

o 
o 

<D 
> 

Oh 

<D 
O 
O 



10 
10 
10 
10 
10 
10 
10 



-1 

-2 
-3 
-4 
-5 
-6 
-7 



o 
o 

C/3 
C/3 

O 

o 

<D 
-I— > 

Oh 

<D 
O 
O 



0.4 
0.3 
0.2 
0.1 








^bcain 


40 MeV 










/ \ V 








































1 1 1 1 1 1 1 1 1 1 1 1 1 







20 



40 



60 80 100 

My (MeV) 



My = 20 MeV 




My = 10 MeV 



M^= 5 MeV 



20 40 60 80 



100 120 140 



"'beam 



(MeV) 



Figure 5: signal cross-sections. TOP: accepted cross-section as a function of Mu- The 
solid line shows the result with a fixed coupling constant, and the dashed line shows the 
result when the constant varies with mass according to Eq. (Qp. Here, '^accepted" refers to 
a limited angular range for the scattered electron: 89.4° < 9 < 90.6°. BOTTOM: cross- 
section as a function of the beam energy, for three values of Mu, as indicated. The curves 
for Mjj = 5 MeV and 20 MeV have been multiplied by factors of 20 and 0.2, respectively. 
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Figure 6: differential signal cross-sections as a function of angles. TOP: da/dQ plotted as 
a function of cos 6. The elastic scattering cross-section has been scaled down by a factor 
of 10"^''. BOTTOM: da/dOp showing the broad distribution expected for the signal. The 
vertical dashed lines show the very narrow range expected for elastic scattering even allowing 
for radiative processes. 
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Figure 7: signal shapes compared to a Gaussian centered at E' = 38.4 MeV with an r.m.s. 
of Q% (see text). Two examples of phase-space distributions are shown, for = 2 MeV 
and 9 MeV (see Sec. \4.3.1\ ). One example of a CompHEP matrix- element calculation is 
shown, for M^ = 2 MeV and Mu = 10 MeV (see Sec. \13^. 



5 Proposed Experiments 

The two proposed experiments bring a low-energy electron beam (E'beam = 40 MeV) onto a 
fixed hydrogen target, and record a large number of elastic scattering events. Each design 
exploits kinematic differences as discussed in Section H] above. The first design is simple and 
is meant to demonstrate the concept, while the second is somewhat more ambitious, though 
easily within the capabilities of existing techniques. 



5.1 Experimental Design I 

This design exploits the correlations between the scattering angle and outgoing electron 
energy, and the proton angle and velocity. Fixing the electron scattering angle in a narrow 
range 89.4° < 6 < 90.6°, the other three quantities should be centered in narrow ranges 
around certain values, namely, scattered electron energy E' = 38.4 MeV, scattered proton 
momentum = 55.4 MeV, and angle 6p = 43.8°. We imagine that the electron is detected in 
a high-resolution electromagnetic calorimeter placed at ^ = 90° and covering 27r in azimuth, 
and the proton is detected in a Silicon strip detector placed at 6p = 43.8°. Elastic scattering 
would give a coincidence in the calorimeter and Silicon detectors. The calorimeter would 
measure the energy of the electron, and the silicon would measure the kinetic energy of the 
proton due to its very high dE/dX, since v ~ 10~^c for the proton. 
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A redundant measurement of the proton velocity would come from measuring its time- 
of-flight across a fixed distance. We present this time-of-flight as a function of the proton 
momentum in Fig. [HI where the dot indicates the expected value = 55.4 MeV for elastic 
scattering. The resolution on the proton momentum worsens as the momentum increases. 
Nonetheless, for a distance d = 2 m, and for a time-of-flight resolution of 2 ns, the mea- 
surement uncertainty on would be only 1 MeV, which is more than adequate for these 
purposes. 

The energy resolution of the calorimeter is of paramount importance. Recent examples of 
high-resolution electromagnetic calorimetry based on high-purity single crystals suggest that 
a resolution of about 6% could be achieved for E' = 38.4 MeV [41]. Fig. [7] compares a purely 
Gaussian peak centered on this value for £" with an r.m.s. of 6% to three signal distributions 
generated according to the phase space model, and the matrix-element calculation. The 
Gaussian curve is normalized to unity and the three signal shapes are normalized to 10"^*^. 

A clear separation between signal and background from the elastic peak is evident, and 
in this simple picture it is easy in principle to isolate a signal by requiring E' < 20 MeV, for 
example. This cut would retain more than 75% of the signal, while the Gaussian is effectively 
removed. 

The signal from elastic scattering will not correspond exactly to a Gaussian, however, 
due to resolution effects in the calorimeter and due to radiative processes mentioned in 
Section m It is beyond the scope of this paper to try to estimate the exact shape of the elastic 
scattering peak. One should note that many of the radiated photons will be emitted along 
the direction of the outgoing electron, and hence will strike the same calorimeter element 
that the scattered electron does. This element will therefore register the energy radiated 
by the electron together with the photon, providing a kind of automatic correction for the 
radiation. If photons are emitted by the incoming electron resulting in an effective radiative 
energy loss for the interaction, the energy of the outgoing electron and proton decrease 
proportionally, as shown in the top plot of Fig. [91 However, the angle of the outgoing proton 
hardly changes at all, as shown in the bottom plot of Fig. [9] Consequently, a robust signal 
for elastic scattering is available. Radiation emitted at wide angles is extremely rare, perhaps 
on the order of 10~^ [37], and in principle could be tagged. 

Scattering from residual gas molecules which have migrated far from the nominal inter- 
action point might mimic the kinematics of the signal, and thereby constitute a background 
that is not easy to reject on the basis of our kinematic cuts. An assessment of the level of 
this background requires a detailed design and perhaps some level of proto-typing, which is 
beyond the scope of this paper. 

It is difficult to estimate the power of this apparatus to reject elastic scattering. Nonethe- 
less, it may be plausible that our signal could be isolated from elastic scattering. The sig- 
nature would be an outgoing electron at low E' values (well below 25 MeV) and either the 
absence of a proton, or a slow proton exiting at an angle significantly different from 44° (see 
bottom plot of Fig. [5]), and with an azimuthal angle with respect to the electron different 
from 180°. QED radiative processes would not give a peak at low E', rather, they would 
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Figure 8: time- of-flight for the outgoing proton, when it traverses a distance of two meters. 
The dot marks the expected momentum for elastic scattering in experimental design I. 



produce a tail to the Gaussian peak shown in Fig. [71 The variation of any signal with -Ebcam 
and 6 is quite different from the variation of the background (see Fig. [5] (BOTTOM) and 
Fig. E] (TOP)). 

As a toy example, we imagine an incoming beam with E = 40 MeV and a cur- 
rent of 10 mA. It impinges on a gas-jet hydrogen target with a density of p = 2 x 
10^^ atoms/cm^ [42]. Our calorimeter consists of one ring of elements arrayed at 6' = 90°, as 
shown in Fig. [TUl The calorimeter elements are single high-purity PbW04 crystals shaped 
as regular polyhedra which project toward the target. The front face might be 2 cm x 2 cm, 
and the elements might be placed 1 m from the target. We assume a Gaussian resolution 
of 6%, or 2.4 MeV at 38 MeV [41]. Additional plastic scintillator would be used to indicate 
when a shower has leaked laterally outside the crystals, and when a wide-angle photon has 
been emitted. The proton is detected with a silicon strip detector placed at the angle ex- 
pected for elastic scattering. Additional detectors could be placed at other angles to confirm 
a proton scattered in the signal process. The experimental region would be filled with a 
low-Z material, such as Helium gas, to minimize multiple scattering and energy loss of the 
proton. 

We integrated Eq. [2] numerically to obtain an accepted cross-section for this simple ap- 
paratus, and obtained 0.8 /ib. With the stated beam current and target density, the rate of 
elastic scatters is 2 x 10^ events/s. With two months' good running time (10 hours per day 
for fifty days), one would accumulate on the order of 10^^ elastic scattering events. 

An accurate estimate of the sensitivity of this experiment to our signal process would 
require detailed simulations of the apparatus, radiative and other energy loss processes. 
We cannot provide such an estimate at this time. However, a rough order-of-magnitude 




Figure 9: the impact of radiative energy loss due to photons emitted by the incoming elec- 
tron on the kinematics of elastic scattering. TOP: the momentum of the outgoing particles. 
BOTTOM: the variation of the proton scattering angle 9p. 
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estimate can be made as follows. Fig. [7] shows a very clear kinematic distinction based 
on the measured scattered electron energy. A simple cut such as E' < 20 MeV amounts 
to 8a below the elastic peak. For an ideal calorimeter, none of the 10^^ elastic scatters 
would survive this, while > 75% of the signal would survive, depending on M^. In practice, 
this simple cut on E' would be reinforced by the powerful cuts on the outgoing proton. A 
positive signal of 15 events would correspond to a signal cross-section that is a about factor 
of 2 X 10"^'^ smaller than the elastic scattering cross-section, or about 0.2 fb. The absence of 
any signal would correspond to an upper limit of about 3 x 10~^ pb, for our default choice 
of coupling constants and masses. We can estimate how such an upper limit would vary as 
a function of Mj/, and the result is shown as the hook-shaped curve in Fig. [51 Apparently, 
this simple apparatus has a good chance to observe a signal or rule out this model. 

If a signal were observed in the energy spectrum of the scattered electron, then one could 
expand the coverage of the Silicon detector so that the recoil proton would be well measured 
for signal events. Then the kinematics of the initial state and the outgoing electron and 
proton would be known on an event-by-event basis, and it would be possible to measure 
the mass of the invisible final state by computing the recoil mass. Specifically, with our 
definition of the missing four-momentum, W^j, = — p'^ 
the two dark-matter particles would be M, 



XX 



— P^, the invariant mass of 
Wf.W^'. There will be a threshold of 2M^ 



for JW^W^^, and if the f/-boson is produced on mass shell, then there will be a peak at 



WfjW^ = Mu. A detailed simulation is needed in order to understand the resolution of 
the recoil mass. 
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5.2 Experimental Design II 

This design augments the calorimetry of the first design with an open-geometry tracking 
device and a magnetic field. Precise measurements of the momenta and direction of the 
outgoing electron and proton would provide powerful handles against background processes. 
As pointed out in Section |U there are angular variables which can be exploited to identify a 
signal. Furthermore, such measurements serve to verify the electron energy as measured by 
the calorimeter and thereby reduce radiative backgrounds. 

We envision a cylindrical geometry coaxial with the beam - see Fig. [HI The region 
around the target would be instrumented with two layers of precision silicon detectors which 
would provide the position and angles of the outgoing electron and proton. For the sake of 
discussion, we take the radii of these two cylinders to be 1.5 cm and 2.5 cm. It is likely that 
the inner cylinder could be reduced to a few mm, depending on the properties of the beam. 
Directly outside the second cylinder we generate a toroidal magnetic field. A small number of 
conducting rods will run parallel to the beam and carry a modest current. A circumferential 
magnetic field is generated with a strength on the order of hundreds of Gauss. This field 
bends the electron and proton but they remain co-planar with the beam, which is important 
for triggering and data analysis. At a radius of 1 m, we install a cathode strip or multi- 
wire proportional chamber to measure the 2;-coordinates of the electron and proton. This 
measurement together with the initial trajectory measured by the silicon detectors gives us 
the momenta of the two particles. Another set of conducting rods carries a current which 
terminates the magnetic field. Finally, an extended version of the calorimeter measures the 
energy of the scattered electron, independently of the tracking system, as well as any photons 
radiated away from the beam. 

All of the proposed detector elements are ordinary, and many examples exist in working 
detectors. The scale of the apparatus is modest, as is the number of instrumented channels. 
The event rate would be rather low compared to modern high-energy physics experiments, 
so read-out electronics would not need to be especially fast. 

5.2.1 Analysis Plan 

The measurements zo,zi, Z2 of the electron position at the surfaces of the three nested cylin- 
ders allow an accurate identification of a dark matter production event. The accuracy of 
this analysis plan depends significantly on the properties of the silicon detectors used at 
the surfaces of the inner two cylinders; in particular, there is an absolute "worst-case" error 
Az in the measurement of position, on the order of 100 /xm. In the following discussion 
we suppose that the incoming electron beam is traveling along the 2;-axis in the positive 
direction. For a given event, we compute the measured scattering angle ^meas from zo,zi, 
and the distance Ri — Rq between the two inner cylinders. We then simulate the path of the 
electron assuming the event was an elastic collision, obtaining the outgoing electron energy 
from the kinematics of the reaction, and find the expected position Z2 at the outer detector 
wires. At this point all events with Z2 < z'^ are rejected. 
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Figure 11: experimental design II 



If the event were in fact an elastic collision, and ^meas > Gtme, then Z2 > z'2 and such an 
event would survive the first cut. To eliminate these remaining elastic events, we consider 
the worst-case scenario when ^mcas is too large; that is, where there are measurement errors 
of +Az at Ro and —Az at -Ri. In this situation, our "worst-case" value of ^tme is computed 
using zq — Az and zi + Az, instead of zq and zi. We then find the outgoing energy for an 
elastic collision with this scattering angle and repeat the simulation to find the new expected 
position Z2. After cutting all events with Z2 < Z2, we are guaranteed that no elastic events 
have survived. For 70° < ^truc < 110°, ^bcam = 40 MeV, and = 2 MeV, a simple Monte 
Carlo simulation shows that the majority of dark matter production events survive the cuts. 



5.2.2 Rates and Sensitivity 

The rate of elastic events would scale up from what was discussed in Section 15.11 due to the 
much larger solid angle covered. Also, one would presumably employ a more intense beam 
and/or thicker target. We imagine that on the order of 10^^ elastic scattering events would 
be recorded, giving a sensitivity on the order of 10"*^ pb. More importantly, a large sample 
of signal events could be detected, allowing a real measurement of the cross-section and of 
the mass of the ?7-boson. 
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5.2.3 Threshold Behavior 

Fig. [5] (BOTTOM) shows the excitation curves for our signal process. Our default beam 
energy, E'beam = 40 MeV, is close to the maximum cross-section when Mu = 10 MeV. 
With the larger event samples expected in this more sophisticated experimental design, one 
could measure the cross-section as a function of -Ebeam; and trace out the threshold behavior. 
Fig. indicates that even just three modest measurements at -Ebeam = 40 MeV, 20 MeV and 
80 MeV would allow one to distinguish clearly between the three mass values Mu = 10 MeV, 
20 MeV and 5 MeV. Furthermore, experimental design II would allow a measururement of 
the angular variation of the signal process. As shown in Fig. [6] (TOP), a clear distinction 
with respect to elastic scattering could be made. 

5.2.4 Other Signatures 

Returning to the light dark matter model, we recall that the intermediate particles ([/-boson 
and/or F^-fermions) couple to electrons, and hence the process 

e~p e~pe^e~ 

must occur (see Fig. [T2|) . The invariant mass of the extra e~^e~ pair would peak at the 
f/-boson mass. If the electron coupling /g is not too low, then one might obtain a sample 
of such events, which would be easy to distinguish from any other process. The observation 
and measurement of all final-state particles would allow an incisive study of this model. The 
prediction for the rate of this process, however, is more model dependent, since it depends 
on /g rather than C^fe, and is less constrained by the dark-matter abundance, Eq. (Q. 

5.3 Comparison to e+e -Scattering 

A classic study of invisible final states has been carried out at e~^e~ colliders. The idea is to 
tag the event by the presence of a wide-angle photon radiated from the incoming electron or 
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positron. Recently, Borodatchenkova et al. proposed using this process to detect light dark 
matter [3]. Also, Zhu discussed the possibility of making an observation at BES III [43], and 
claims a sensitivity down to /e ~ 10~^. 

Our proposal benefits in a number of ways. First, the luminosity is much higher due 
to the use of a fixed target. Second, the kinematics are particularly advantageous, since at 
tree-level there are no backgrounds and elastic scattering is highly constrained. Finally, the 
recoil mass can be used to measure the f/-boson mass, in principle. 



6 Summary and Conclusions 

The light dark matter model proposed by Boehm and Fayet can explain both the relic density 
of dark matter and the 511 keV gamma-ray line coming from the center of the galaxy. It 
survives many constraints derived from precision measurements in particle physics, and also 
several astrophysical observations. The key ingredients of this model include a stable light 
neutral scalar, x, and a light neutral vector boson U with a large coupling to x ^"^^ ^ small 
coupling to electrons. 

We propose a simple experiment to produce the fZ-bosons in low-energy electron-proton 
scattering. Basically, the t/-boson would be radiated from the incoming electron, resulting 
an invisible final state which can be identified through dramatic kinematical differences with 
respect to elastic scattering. The high luminosity and easy control of the kinematics should 
allow the experiment to reach a sensitivity which would confirm or definitively rule out this 
model. Two conceptual designs are described, which will be modeled and studied in greater 
detail in the future. 
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